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Abstract 

We investigate the nonlinear holomorphic supersymmetry for quantum-mechanical 
systems on Riemann surfaces subjected to an external magnetic field. The realization 
is shown to be possible only for Riemann surfaces with constant curvature metrics. 
The cases of the sphere and Lobachevski plane are elaborated in detail. The partial 
algebraization of the spectrum of the corresponding Hamiltonians is proved by the 
reduction to one-dimensional quasi-exactly solvable s((2,R) families. It is found that 
these families possess the "duality" transformations, which form a discrete group of 
symmetries of the corresponding ID potentials and partially relate the spectra of dif- 
ferent 2D systems. The algebraic structure of the systems on the sphere and hyperbolic 
plane is explored in the context of the Onsager algebra associated with the nonlinear 
holomorphic supersymmetry. Inspired by this analysis, a general algebraic method for 
obtaining the covariant form of integrals of motion of the quantum systems in external 
fields is proposed. 
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1 Introduction 



In the pioneer paper QXJ] , the usual linear supersymmetry [§] was generalized by employing the 
higher-derivative supercharges. The characteristic property of such a generalization is the 
polynomiality of the corresponding superalgebra in even integrals of motion of the system. 
This makes the polynomial (nonlinear) supersymmetry to be similar to the Yangian and 
finite W-algebras || [|, |5| . The supersymmetry of such a type was found in various physical 
models || [7|, [| [TO, O. This provides a solid background for physical interest in the 
nonlinear supersymmetry. 



The nonlinear holomorphic supersymmetry (n-HSUSY) JL2J, [13], [TJJ |15| is a natural gener- 
alization of the usual linear || and higher- derivative (polynomial) [I], ^, [l(j] supersymmetries. 
Its construction was triggered by the observation of the quantum anomaly problem, which 
appears under attempt to quantize the classical one-dimensional system with nonlinear su- 
persymmetry of arbitrary order n 6 N, n > 1 |. The important result in resolving this 
problem was obtained in Ref. [121, where we showed that the anomaly- free quantization 
is possible only for the peculiar class of the superpotentials, for which the corresponding 
quantum n-HSUSY turns out to be directly related to the quasi-exactly solvable (QES) 
systems [[L6|, [L7|, 



18 , |19[| . Such one-dimensional quantum mechanical systems were also inde- 



pendently discussed within the framework of the so called Type A TV-fold supersymmetry 
PH pl[ |22| , p3| , |2~3| , where, in particular, the equivalence between the nonlinear supersym- 
metry and the sl(2,R) scheme for the ID QES systems was demonstrated |22 |. 

On the other hand, nowadays, the two-dimensional dynamics attracts a considerable 
attention in the context of the both, linear 



25 



and nonlinear 28, 29, 30 



su- 
persymmetries. The particularly interesting 2D system, related to the quantum Hall effect 
f31[ ) [32], is the non-relativistic charged spin-1/2 particle in a stationary magnetic field 
|3~4] , [351 |25| . In Ref. jl3| , we generalized the n-HSUSY to the 2D case represented by the 
system of a charged spin-1/2 particle with gyromagnetic ratio 2n (or, spin-n/2 particle with 
gyromagnetic ratio 2) moving in the plane with magnetic field of a specific form. The most 
important result of such a generalization was the observation of some universal nonlinear 
algebraic relations underlying the n-HSUSY. 

In the recent paper flT5| , the nonlinear relations of Ref. [|13| were identified as the Dolan- 
Grady relations p6| . As a result, the construction of the n-HSUSY was reduced to the 
following universal (representation-independent) algebraic structure. Its supercharges have 
the form of (anti) holomorphic polynomials in the two mutually conjugate operators (gener- 
ating elements), in whose terms all the components of the n-HSUSY construction are built. 
It is these generating elements that obey the Dolan-Grady relations and produce the asso- 
ciated infinite-dimensional Onsager algebra |[7). The knowledge of the universal algebraic 
construction essentially facilitated the search for the central charges of the n-HSUSY, and 
simplified the calculation of its nonlinear superalgebra. Having it, we proposed a generaliza- 
tion of the n-HSUSY in the form of nonlinear pseudo-supersymmetry |T5]] , which is related 
to the PT-symmetric quantum mechanics |38| , |39| , (40| . 

The power of the algebraic formulation of the n-HSUSY [JOJ is that it does not refer to 
the nature of the space-time manifold which can be lattice, continuum or noncommutative 
space. In this sense the n-HSUSY is a direct algebraic generalization of the usual quantum- 
mechanical supersymmetry H. 
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The present paper is devoted to investigation of the quantum mechanical problem of 
realization of the n-HSUSY on Riemann surfaces with an inhomogeneous magnetic field. 



Applying the algebraic construction of Ref. ||15||, we find a necessary condition for such a 



realization, and discuss the spectral problem in the two possible (non-trivial) cases of the 
sphere and Lobachevski plane subjected to an external axial magnetic field of a peculiar form. 
Reducing the n-HSUSY systems to one dimension, we find the new discrete symmetry of 
the ID QES potentials, which induces the "duality" transformations between the associated 
distinct sl(2, M) schemes. The analysis of the axial symmetry of the 2D n-HSUSY systems 
in terms of the intrinsic algebra generated by the covariant derivatives allows us to propose 
(beyond the context of supersymmetry) the general algebraic method for obtaining the co- 
variant form of integrals of motion of quantum systems in a curved space in the presence of 
an external gauge field. 

The paper is organized as follows. In Section 2 we discuss the realization of the n-HSUSY 
on Riemann surfaces with a magnetic field, and find the zero modes corresponding to the 
cases of the sphere and Lobachevski plane. In Section 3 we demonstrate that the spectral 
problem for these two cases can be partially solved by a reduction to one dimension. The 
obtained ID systems admit distinct s[(2,R) representations, that is reflected in existence 
of the discrete symmetry of the potentials. This symmetry is investigated in Section 4. In 
Section 5 we study the algebraic content of the systems on the sphere and hyperbolic plane 
in the context of the associated infinite-dimensional contracted Onsager algebra ||15|| . Section 
6 develops the algebraic approach for finding the covariant form of the integrals of motion 
of the quantum systems in an external field with symmetry. In Section 7 we discuss the 
obtained results and specify some problems to be interesting for further consideration. 

2 n-HSUSY and Riemann surfaces 
2.1 n-HSUSY algebraic structure 



The system with nonlinear holomorphic supersymmetry is described by the Hamiltonian [15| 



! — 77 — 

H n = -{Z, Z}+-[Z, Z]a 3 . (2.1) 

Here 03 is the diagonal Pauli matrix, n e N, and the mutually conjugate operators Z and Z 
obey the nonlinear Dolan-Grady relations 

[Z, [Z, [Z, Z]]]=u 2 [Z, Z] , [Z, [Z, [Z, 2}}} = uj 2 [Z, 2} . (2.2) 

These relations guarantee the existence of the odd integrals of motion (supercharges), Q n 
and Q n , which are defined by the recurrent relations 

Q n =(z 2 -(^) 2 u 2 )Q n _ 2 , Q = v+, Qi = Za + , 

) ' (2-3) 
Q n = [Z 2 - (2fi) V) Q n _ 2 , Qo = <r-, Qi = Za., 

containing a constant parameter to G C, uj = uj*, and a± = \((J\ ± ia 2 ). The supercharges 
Qn, Qn generate a nonlinear superalgebra of order n []15[|. In Section 5 we shall return to the 



3 



discussion of the algebraic construction underlying the nonlinear holomorphic supersymme- 
try. 



2.2 Fixing geometry 

We are going to investigate the problem of realization of the n-HSUSY on a 2D surface 
with a nontrivial Riemann geometry. The metrics of any two-dimensional surface can be 
represented locally in the conformally flat form, 

ds 2 = p(z, z)dzdz = g zz dzdz + g zz dzdz, (2.4) 

where z = x 1 + ix 2 , z = x l — ix 2 are the isothermal coordinates. Given the isothermal 
coordinates on the neighbouring patches, the metrics on the patches are related by conformal 
transformations. Hence, we have 

p - 2 

9zz = -> 9 ZZ = ~, 9zz = 9zz = 0, V9 = P, 

2 p 

where g = detg^. 

Let us choose the operators Z and Z in the form of covariant derivatives, 

Z = V 2 = D + d log p ■ S, Z = -V 2 = —g z2 D, (2.5) 

where the operator S defines the helicity, or the conformal spin of a tensor T^, ST^ = sT^ s \ 
s£Z, while D = d + ieA and D = d + ieA are the U(l) covariant derivatives with A, A to 
be the components of the vector potential corresponding to the case of stationary magnetic 
field. The operator e gives the electric charge of the tensor field. By definition, the derivative 
V z maps the tensor fields of the helicity s into those of the helicity s — 1, while V z transforms 
the tensors of the helicity s into those of the helicity s + 1. The operators ( [2.5D are mutually 
conjugate with respect to the scalar product [f|l|] 

(T^ S \U^) = J d 2 z^(g zS ) s f^U^ (2.6) 

defined on the space of the tensor fields satisfying the appropriate boundary conditions. 
In representation ( |2.5| ) we have 

[Z,Z] = Be-~llS, (2.7) 
where the scalar curvature of the Riemann surface TZ and the magnetic field B are 

K = —ddkip, Be = --[D,D]. (2.8) 
P P 

As a result, the nonlinear Dolan-Grady relations (|2.2|) take the form 

((V 2 Z B) - uo 2 B) e - ~ ((V^) - u 2 n) S - (V Z H)V Z = 0, 
{{V z2 B) - u 2 B) e-\ ({V z2 TZ) - cu 2 K) S + (V z ft) V 2 = 0. 



(2.9) 



4 



The relations (2.S) have to be satisfied as the operator identities. Hence, the coefficients at 
the operators e, S, V z and V 2 have to vanish independently. The coefficients at the covariant 
derivatives vanish only when the scalar curvature is a constant, 

K = const. (2.10) 

With this condition, the coefficients at the operator S result in the equation u) 2 7Z = 0, and, 
hence, there are two possibilities: either uj = 0, or 1Z = 0. The latter case corresponds to 



the plane geometry, and it was analysed in detail in Ref. Therefore, in what follows we 
shall discuss the former case. 

For lu — 0, there arise the contracted Dolan-Grady relations, 

[Z, [Z, [Z, Z]]] = 0, [Z, [Z, [Z, Z]]] = 0. (2.11) 



The nonlinear supersymmetry associated with the contracted Dolan-Grady relations Q2.ll 
is generated by the supercharges of the form 

Q n = Z n a + , Q n = Z n a^, (2.12) 



corresponding to the limit u — > of the general case ( |2.3| ). 

One notes that, in principle, the restriction (|2.10|) can be overcome if to endow the 
Riemann structure of the 2D surface with the torsion. Then, for u ^ 0, the Dolan-Grady 



relations (|2.2| ) will result in a set of nonlinear equations on the scalar curvature and on the 
torsion. However, here we restrict ourselves by the torsion-free case with the relation ( |2.10|) 
to be the necessary condition for existence of the n-HSUSY on a Riemann surface. 



The general uniformization theorem [42j implies that, up to the holomorphic equivalence, 
there are just three distinct simply connected Riemann surfaces: a) the sphere S 2 with the 
conformal factor 

P(?>z) = n | L -x 2 > ( 2 - 13 ) 
(1 + p z zz) 2 

where z G C U {oo}; b) the complex plane C with the flat metric, p = 1; c) the Lobachevski 
(hyperbolic) plane with the metric ( |2.4| ) defined on the disk \z\ < fi~ x by 



In Eqs. (p.l3|) , (|2.14|) , the positive parameter (3 defines the value of the scalar curvature, 



1Z = ±8/5 2 , where the plus corresponds to the sphere, while the minus does to the hyperbolic 
plane. In the case of the sphere the parameter (3 is related to its radius as = 2R. 

In what follows we will mainly analyse the case of the sphere, while all the corresponding 
results for the Lobachevski plane can formally be obtained by the change (3 — > i/3. At the 
same time, the construction of the n-HSUSY on the plane [TB|] can be reproduced from the 
cases corresponding to the surfaces ( |2.13| ) and ( 2.14j) via the appropriate limit procedure 
(3^0. 
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2.3 Fixing magnetic field 

We have not analysed yet the conditions corresponding to the vanishing of the coefficients 
at e. They produce the equations on the external magnetic field. In correspondence with 
Eq. (|2.8|) , the magnetic field is given locally by the relation 

Bdv = dA, 

where dv = dz A dz is the area (volume) element, and A = Adz + Adz. In cor- 

respondence with ( |2.9|) , for uo = we have the following equations on the magnetic field 
B = -ig zS (dA -OA): 



For the both cases ( |2.13|) and (|2.14j ), their solution can be represented in the form 



B(z,z) = c 1 ^Jp{z,z) + c , (2.15) 
with Co, C\ G K. In the case of the sphere the (normalized for 27r) total magnetic flux, 

oo 

27r$ = J Bdv = 2tt J B(r)p{r)rdr = ^g|_t^l ) ( 2 .16) 

s 2 o 

is quantized, $6Z. Here we have used the polar representation of the complex coordinates, 

z = re lip , z = re~' lip . 



In the case of the hyperbolic plane, the magnetic flux is divergent, but the relation ( |2.16| ) 



(with the substitution (3 2 — > — (3 2 ) still can be used in the form $ = — (2c + Ci)/(4/? 2 ) as a 
mere combination of the parameters which can acquire any real value. 
In what follows, we fix the charge of the system to be e = +1. 

2.4 Zero modes 



Let us find the zero mode subspace of the n-HSUSY realized on the Riemann surfaces ( p,13| ) 



fl2.14| ) with magnetic field ( |2.15|) . To this end we note that the space of the states of the 



system ( |2.1| ) is a direct sum of the two Hilbert spaces, JF_©JF + , formed by the states G T- 
and G being the eigenstates of 03, a^ 1 " = ^f^/ T . We will refer to the spaces T- 
and respectively, as to the " bosonic" and "fermionic" subspaces. We also assume that 
the "bosonic" wave functions are the tensors of helicity s. Then, in accordance with the 



structure of the supercharges ( |2.12|) , the "fermionic" wave functions will be the tensors of 
the helicity s — n. 

The Hamiltonian acting on the "bosonic" wave functions can be rewritten in the complex 
coordinates as 

#W= -j-({D,D} + 2sd\ogp-D)-^B-s(n + l)p 2 . (2.17) 
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The upper index in the parentheses denotes the helicity of the tensor that the Hamiltonian 
acts on. In these notations the "fermionic" Hamiltonian is ifi s ~ n ' ) . 
In the gauge 

A = -i Z - (c lP + 4/5 2 $Vp) , A = z| (cip + 4/3 2 $ v /p) , (2.18) 

corresponding to the magnetic field ( [2.15| ), ( |2.16|) , the states of the form 

n-1 

^-\z,z) = P i- s+s ^eW^J2vV(z) zk , (2-19) 



fc=0 
n-1 



k=0 



are the zero modes of the supercharges Q n and Q n , respectively. In the case of the sphere, 
the condition of normalizability fixes the form of the functions tp^\ 

2(n-s $ -l)-fc 2(s*-l)-fc 

^" ) (^)= £ <4^, 4%) = £ ^> ( 2 - 2 °) 

q=0 9=0 

where we have introduced the notation 

s<& = s __$. (2.21) 

The physical sense of the parameter s$ will be found below under analysis of the contracted 
Onsager algebra associated with the system. The corresponding number of the normalizable 
zero modes in the "bosonic" and "fermionic" sectors is given by 

n b = (3n - 4s$ - 1) , Np = — (n - 4s$ + 1) . 

Here we assume that Ns,Np £ N; for other values of the parameters n, s and $ the 
corresponding number of zero modes is implied to be equal to zero. One can verify that in 
the spheric case iV^ + iVjr > for any n £ N, and, therefore, the nonlinear supersymmetry 
of the corresponding system is always unbroken. 

In the case of the hyperbolic plane, the normalizability is managed by the sign of the 
parameter ci, and the corresponding space of zero modes is infinite-dimensional. Moreover 
on the hyperbolic plane, unlike the sphere case, the supercharges Q n and Q n can not have 
zero modes simultaneously. 

So, the nonlinear holomorphic supersymmetry of quantum-mechanical systems on 2D 
surfaces leads to the strong restrictions on configurations of the external geometry and mag- 
netic field. Moreover, the systems with the n-HSUSY on the sphere are quasi-exactly solvable 
since the space of zero modes of the supercharges is finite-dimensional in this case. In the 
next section we will demonstrate that the same is also true for the systems on the hyperplane. 
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3 Reduction to s[(2,R) quasi-exactly solvable families 



In this section we consider the spectral problem for the "bosonic" Hamiltonian ( [2.171 ) an d 
show that it can be reduced to a one- dimensional problem corresponding to s[(2,M) quasi- 
exactly solvable families of the systems. The spectral problem in the "fermionic" sector can 
be analysed in a similar way. 

First, we note that the functions of the form z m f(zz), m G Z, constitute a subspace 
invariant with respect to the action of the Hamiltonian (|2 . 1 7| ) in the gauge (|2.18|) . The 
parameter m, as will be shown, is the quantum number corresponding to the integral ( |5.6| ) 
associated with the axial symmetry of the system. Departing from the form of zero modes 
( p,19| ), we search for the eigenfunctions of the form 



V>(z, z) = z m {zz) a {p{zz)) b e c ^ ir) <p(zz). (3.1) 
Now, we introduce the variable 

with the domain < y < 1. The anzatz ( |3.1| ) reduces the two-dimensional Hamiltonian 
( p. 17 ) to the one-dimensional operator 

H s i(2,m) = 2/3 2 (Tq - T_T ) + a + T + + a T + a„T_ + const (3.2) 

if the parameters a, b and c are chosen in the following way: 

a = m, o = —s<f, m, c = u—^. (3.3) 

2 2 2 4' P 8/5 2 V ; 

The parameters /z, v and a acquire the values ±1, while the operators 

T + = {N -l )y -y*-, T = y--— r. = _ (3.4) 

generate the sl(2, M) algebra with the Casimir operator fixed by the parameter N. When 
N G N, the realization (|3.4| ) gives an irreducible N- dimensional representation on the space 
of polynomials of the degree up to iV — 1. In this case the operator ( |3.2p corresponds to some 
QES system ITH, [1] 



The corresponding parameters of the sl(2, M.) scheme, N, a± and «o, read as 

-m — (/i + ^)s$, 



N 


= /in H 




Cl 


a + 




a 






= f? (l- 



3u + a , 
fin H — m + [(/, — 3z/)s$ I . 



So, we have formally 2 3 = 8 solutions parametrised by the triple (/i, z/, a). Of course, the 
solutions (|3.5|) realize a true quasi-exactly solvable system only if the condition iV G N is 



S 



satisfied. For example, the solution with the triple (—1, +1, +1) does not correspond to any 
QES system since in this case N = —n, while throughout the paper we assume n e N. 
But as it will be shown further, this solution together with the other solutions reflects some 
discrete symmetry of the resulting one-dimensional system. We also note that under the 
change (fj,,u,cr) — ► (—fjL,—u,—a) in (|3.5| ), the parameter N changes its sign. This means 
that for any choice of the parameters n, m and s$, at most four solutions from Q3.5Q can 
realize the proper finite-dimensional s[(2, R) representations of the QES system. 

From the 2D viewpoint, the conditions of regularity at zero and of finiteness at infinity^ 
for the functions (|3.1|), ( |3.3j ) lead, correspondingly, to the constraints 

am > 0, urn < 2(vs$ — s). (3.6) 

The condition of normalizability of the wave functions results, in turn, in the inequality 

z/(2s$ -m) + 1 > 0. (3.7) 

In all the relations ( |3.6|) , (|3.7| ) we imply that the parameters m, s and <3> accept integer values 
only. 

Changing the variable, 

y = cos 2 (3x, 

with < x < and following the standard procedure of the s((2,R) scheme for QES 
systems [If], 17, 18], one can transform the operator (|3.2j ) to the one- dimensional Hamiltonian 
of the standard form, 

All the solutions Q3.5p lead to the very potential 



V(x) = - ^cos 4 (3x + | f 2 (a, - n) + ^ j cos 2 (3x 

+ tL. (Am 2 - l) cot 2 fix + ^- (4 (2s $ — m) — l) tan 2 (3x. 



(3.J 



This means that the potential (|3.8[) possesses some discrete symmetry, which relates the 
solutions ( p.5| ). We discuss this symmetry in the next section. 

The potential ( |3.8|) has the following behaviour near the point x = 0: 

f , 4m 2 — 1 „ . 9 . 

V(x) = const H h C(a; ). 

8x^ 

Therefore, in the case m = ±|, it is formally possible to extend the domain of definition of 
x to the symmetric interval — ^ < x < ^ (or to R in the case of the hyperbolic plane, see 
below) . 

The case of the plane corresponds to the limit /3 — ^ 0. Rescaling and shifting the param- 
eters, 

. Ci Ci 

c -c + ^, c ^~Jv 



1 Since we consider the sphere, it is necessary to require the regularity of all the wave functions for \z\ — > oo. 
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in this limit we recover from Eq. ( j2.15| ) the magnetic field of the form [|13| 

B{r) = C\r 2 + c , r 2 = zz, 

and find that the described reduction results in the one- dimensional QES potential 

V {x) = ||x 6 + ^x 4 + I (eg - 2 Cl (2n - m)) x 2 + + const, (3.9) 

62 o o oX 

defined on the half-line < x < oo. Hence, the proper limit procedure reproduces correctly 



the results derived earlier in Ref. 13 



The formulas corresponding to the case of the hyperbolic plane are obtained by the formal 
substitution (3 — ► i(3 with the change of the domains of the variables y and x to 1 < y < oo 
and < x < oo. Under such transformations, the trigonometric potential (|3.8| ) is converted 
into the hyperbolic QES potential 



V(x) = ^^ cosh4 ( 3x + j[ 2 - cosh 2 fix 

3 2 B 2 
+ — (Am 2 - l) coth 2 fix + — (4 (2^$ - m) 2 - l) tanh 2 (3x. 



(3.10) 



The condition of regularity at zero is the same as for the sphere, am > 0, while the condition 
of finiteness at infinity has to be discarded. Besides, the normalizability of the functions ( |3.1| ), 
( [3. 3D is changed for /iCi > instead of (|3.7| ). Since the magnetic flux diverges in this case, 



the parameter combination $, specified at the very end of Section 2.3, can acquire any real 
value. Therefore, in general, only the solutions (+1, —1, +1), (+1, —1, —1) and (—1, +1, —1) 
can produce the proper finite-dimensional s[(2,lR) representations, while other values of the 
triple (/i, u, a) serve only for integer values of the parameter <3>. 



4 "Duality" transformations 

In the previous section we have seen that the reduction of the two-dimensional Hamiltonian, 
associated with the n-HSUSY on the sphere or hyperbolic plane with magnetic field ( p,15|) , 
results in the one-dimensional quantum system with the potential ( |3.8| ) or ( |3.10| ), respec- 
tively. The resulting ID Hamiltonians can be transformed to the QES operator ( |3.2|) with 
appearance of the eight different families of the s[(2,M) schemes. Here we investigate the 
question what these different families of the s[(2, M.) schemes mean from the point of view of 
the corresponding one-dimensional quantum mechanical QES system. 

First, we note that the potential (|3.8|) is effectively four-parametric since the parameter 
(3 can be absorbed by rescaling x — > (3~ 1 x 1 c\ — > /3 2 ci, H — > [3~ 2 H, while s and $ enter only 
in the combination ([2.21|) . The parameters n, 2s$ and m take integer values, and C\ is real. 
These parameters enter into the potential (|3.8| ) in quadratic combinations. This motivates 
us to search for the linear transformations of the parameters, which leave the potential to 
be invariant. It is convenient to represent them in the matrix form 

Y' = GY, (4.1) 
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with Y T = (n, 2s$, m, c\). Let us introduce the following matrices: 

/-I 





1 

1 

1 








-v 



G-2 



(\ 










G, 



-1 1 
-1 2 
1 

\ -lj \oooi/ 

One can verify that the potential ( |3.8|) stands invariant under the transformations 
generated by the set of the mutually commuting matrices 



(I 





-1 


o\ 





1 


-2 











-1 













V 



(4.2) 



G — 1, -1, G u -Gi 



1,2,3, 



satisfying the relations 



G1G2G3 — — I, 



G\ 



G 2 



G\ 



1. 



(4.3) 



(4.4) 



Here, the eight transformations produced by the matrices G on the space of the parameters 
Y correspond to eight transformations in terms of the parameters (fi, u, cr), which connect 
different solutions ( |3.5|) at the level of the sl(2, M.) schemes. In particular, the transformations 
generated by the G±, G2 and G3 correspond, respectively, to the transformations (/i, z/, a) — > 
(— /i, v, a), (fi, —v, a) and (/i, z/, —a). From this point of view, the transformations ( flTl) -( [Q|) 
give some irreducible representation of the Abelian discrete group of reflections in three axes. 

Another nontrivial example of such discrete transformations is given by the QES potential 
( ^.9[ ) related to the s[(2, M)-operator 

H QES = - 2T T_ + /3 + T + + f3 T + /3_T„ , (4.5) 

where the sl(2, R)-generators are given by Eq. ( |3.4|) with y = x 2 fi~6| , 43| . The complete 
set of the discrete transformations ( f4.1| ) with Y T = (n, m, cq, c{) is given in this case by the 
matrices 1, —1, G, —G, where 



G 



(I 






-1 
-1 






0\ 



1 
1/ 



In comparison with ( |3.8| ), the case ( pT9| ) with n,m e Z, Cq,Ci G K, is more simple, and, 
hence, more illustrative. Let us discuss here the corresponding eigenfunctions, which can be 
found algebraically. One can represent them formally as 



ij(x) = P N ^ (x 2 ) x$ +vm exp (-^c x x 4 - ^c x 2 



with N 



-m, 



(4.6) 



and fx, v acquiring the values ±1. For the functions ( 4.6|) to be normalized, one has to impose 
some restrictions on the parameters. There are four distinct cases, which can be marked by 
the pair (//, v): 



+1,-1) 


N = 


nGff, 


ci > 0, 


rn 


G Z_; 


+1,+1) 


N = 


n — m G N, 


ci > 0, 


rn 


G Z+; 


-1,+1) 


N = 


-n G N, 


ci < 0, 


m 


G Z+; 


-1,-1) 


N = 


m — n G N, 


ci < 0, 


m 


G Z_. 
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The restrictions on c\ arise from the requirement of the normalizability, while those on m do 
from the requirement of the vanishing of the functions at x = 0, since we consider the spectral 
problem on the half-line, < x < oo. In fact, the cases (+1,-1), (— 1,+1) correspond to 
the conventional st(2,R) representation of the potential ( |3.9|) (see, e.g., Refs. |17|, |43[]), while 
those (+1,+1), (—1,-1) were not discussed earlier. All the four cases are related by the 
discrete transformations as it is reflected on the diagram: 



-i,+D <- 



G 



-I 

1,-1) 



In principle, one can treat the system fl3.9| ) from the more general point of view by as- 
suming that all the parameters n, m, Cq and C\ are real. In this case the Hamiltonian can also 
be reduced to the form (|4.5|) , but if there are no additional restrictions on the parameters, 
the sl(2, R) generators realize an infinite-dimensional representation of the algebra. As a 
consequence, the ID system with the potential ( |3.9|) is not quasi-exactly solvable. Now, if 
any of the restrictions ±n 6 Nor ±(n — m) G N is imposed, the ID system under consid- 
eration becomes to be quasi-exactly solvable since with such a restriction the corresponding 
representation of the sl(2,R) is finite-dimensional. For example, when m = |, the domain 
of definition of the potential ( |3.9|) can be extended to the whole real axis. Let c\ > 0, then 
the pair (+1, —1) with N = n £ N gives even eigenfunctions ( |4.(j| ), while the pair (+1, +1) 
with N 

two distinct forms of the QES x 6 -potential [[17|] 

Having the two examples, one can suppose that the existence of the discrete symmetry 
transformations is general for all the one-dimensional QES sl(2, R) systems. Let us give some 
simple arguments in favour of this conjecture. 

Changing appropriately the variable and realizing the similarity transformation, one can 
represent any QES sl(2,R) Hamiltonian in the canonical form [[16], [T7|, 0, [19| 



n — j G N provides the odd eigenfunctions. These two cases correspond to the 



H, 



QES 



N — 1 ., . . A d N , . 

—n(v)-p,(v)) Ty + ^[n(v)- 



N - 1 



P'liv) , (4-7) 



where Pk{y) is a polynomial of the fcth degree. 
Schrodinger equation is given by 



The potential of the corresponding 



V(x) 



N(N + 2) f3P, 
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)/2 



4 P 4 



- P" 



+ 



+ 1 



P' 

2 p' - p 2 t± 
2 Pa 



+ 



1_P[ 

4:P d 



+ const, 



(4.8) 



where the right hand side is evaluated at y 
elliptic integral 

m = 



f 1 (x) with / 1 being the inverse function of 



dy 



'PM 



The canonical form (|4 

which allows us to fix essentially the polynomial P4 



is not unique because of the existence of a "residual" symmetry, 

Wh. Therefore, one can consider that 



12 



coefficients of the polynomial P2 and the discrete parameter iV span the whole parametric 
space of the system. From ( |4.8|) it follows that the parameters enter into the potential in 
quadratic combinations (see the explicit set of QES potentials, e.g., in Ref. (43]). Therefore, 
appealing to the cases (|3.8| ) and ( p.9|) , it seems plausible that the transformations of the form 

Of course, in general, for the transformed 



( |4. 1| ) should also exist in the general case 
parameter N to remain integer, one has to treat some parameters as to be also integer. 

Let us stress that if one treats a ID QES system in the framework of the corresponding 
2D supersymmetric system, then the discrete transformations can intertwine the distinct 
sectors of the spectrum of the same 2D system, or the spectra of different 2D systems. This 
property as well as the relation G 2 = 1 allows us to treat ( |4.1| ) as some kind of "duality" 
transformations. For example, one can verify that the zero modes ( [2. 19 ), ( 2.20| ) are spanned 
by the solutions (|3.3| ) with the triples (+1, —1, +1) and (+1, —1, —1), which are related by 
the "duality" transformation generated by the G3 from Eq. ( |4.2| ). 

In conclusion of this section we note that in the context of QES systems, some similar 
duality transformations were discussed in Refs. 



44, 45, 46 



However, the "duality" transfor- 

(4311261 



mations discussed here are essentially different from the transformations of Refs. 
In the latter case the duality connects different ID QES systems, while in the present case 
it relates different parts of the spectrum of one and the same ID QES system. These parts 
correspond to different s[(2,M) schemes of the given potential. 



5 The contracted Onsager algebra 



In this section, following Ref. P3| , we analyse the algebraic structure underlying the 2D 
spherical and hyperbolic systems with the nonlinear holomorphic super symmetry. 

The operators Zq = Z and Zq = Z together with the contracted Dolan- Grady relations 
( P- 11| ) recursively generate the infinite-dimensional contracted Onsager algebra: 

[Zk, Zi\ = -Bfe+z+i, [Z k , B t ] = Z k+ i, [B k , Z[\ = Z k+ i, ^ 

[Z h ,Z l \=Q, [Z k ,Z t ]=0, [B k ,Bi\ = 0, 

where k, I G Z + and B = is implied. In general, the algebra ( |5.1|) admits the infinite set 



of the commuting charges [15 



1 1 

A = 5 E tt Vi> Z i-v) - BpBi-p) + ? S,(7 S) (5.2) 



2 
P =i 



which includes the Hamiltonian (|2.1| ), J\ = 2H n . 

The operators Z, Z in the representation ( j2.5|) with the conformal factor (|2 .13 ) corre- 
sponding to the sphere generate the following finite-dimensional algebra which we call the 

intrinsic: 

[Z, Z] = G- A(3 2 S$, 
[Z, G] = -4(3 2 V, [Z, G\ = A(3 2 V, [Z, V] = G, [Z, V] = -G, 
[5$, Z) = -Z, [5$, Z) = Z, [5$, V] = -V, [5$, V] = V, (5.3) 
[Z, V] = [V, Z] = [V, V] = [S*, G] = [G, V] = [G, V] = 0, 
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where 

= § - G — B(zz) - 2/3 2 $. (5.4) 

The operator G satisfies the obvious relation jj s2 Gdv = 0, whereas the multiplicative 
operators T> and T> are, in correspondence with Eq. (|5.3|) , the tensor fields of helicity — 1 



and +1, respectively, and have the coordinate representation V = ^cizp(zz), V = \c\z. 
The operators T> and T> are mutually adjoint with respect to the scalar product Q2.6| ). The 
intrinsic algebra ( |5.3| ) and all the corresponding relations for the hyperbolic plane ( 2.14 ) can 
be reproduced by the formal change (3 — > i(3. 

The intrinsic algebra (|5.3|) has the following two Casimir operators: 



C = G 2 + 8p 2 VV, W = VZ + ZV-VV-G{S$-1). (5.5) 

In the given representation, we have C = \c\ ■ J, while the operator W in the gauge ( |2.18|) 
has the form 

W = X -c x (zd -zd- 5$) . (5.6) 



From the coordinate representation ( |5.6| ), one can conclude that the operators D and T> 
represent the components of the Killing vector associated with the axial symmetry, and the 
Casimir operator W is proportional to the generator of this symmetry. It is interesting to note 
that in terms of the algebra ( |5.3| ), the symmetry generator itself has the form J = W/VC, 
i.e. it defines "spin" of representations of the algebra (see below), and the eigenvalue fl2.21p 
of the operator 5$ has a sense of "effective helicity" shifted by the quantized magnetic flux. 
Therefore, in the case of the odd values of the magnetic flux, we have here an example of 
the "boson-fermion" transmutation. We shall return to the discussion of this aspect of the 
system below. 

Rescaling appropriately the generators, it is always possible to reduce the constant 4/3 2 
to e = ±1, where '+' corresponds to the sphere, while '— ' does to the hyperbolic plane. 
Then, one can redefine the generators of the algebra (|5.3|) as 



E„ = Z-^D, E + = Z-^V, E = eSz, (5.7) 

T_ = V, T+ = V, T = eG. (5.8) 

From ( |5.3| ) it follows that the operators ( |5.8| ) span the Abelian algebra of translations, t(3), 
while the operators (|5.7|) together with (|5.8|) obey the commutation relations: 



[E + , E_) = eE , [E , E ± ) = ±E ± , [E , T ± ] = ±T±, 
[E ± , T ] = ±T ± , [E ± , T T ] = T eT , [E ± , T±] = [E , T ] = 0. 

Therefore, one can conclude that in the case of the sphere the intrinsic algebra (|5.3| ) is the 
algebra of the 3D Euclidean group of motions, iso(3) = t(3) Q) s so (3), while in the case of 
the hyperbolic plane it is the Poincare algebra, iso(2, 1) = t(3) © s so (2, 1). In this context, 
the first Casimir operator in (|5.5|) is the "squared mass" operator, C = Tq + 2eT_T + , while 
the second one is the Pauli-Lubanski "pseudo-scalar", W = eT (l — E ) + T + E_ + TJ + T_. 
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In terms of the algebra (|5.3| ), the generators of the contracted Onsager algebra (|5.1| ) read 

as 



Z k = {-Af3 2 ) k {Z + kV), Z k = (-4/3 2 ) fe (Z + kV), 

B k = {-Af3 2 ) k ~\kG-Af3 2 S<,), 

Hence, we have obtained the nontrivial realization of the contracted Onsager algebra in the 
sense that though the infinite set of generators ( |5.9| ) is represented linearly in terms of the 
finite number of the intrinsic algebra generators, nevertheless, Z k , Z k and B\ do not turn 
into zero for k > 1, I > 2. In this context one notes that in the systems with nonlinear 



holomorphic supersymmetry investigated earlier |fL2] , [13| , only the non-contracted Onsager 
algebra case was realized nontrivially ]To| . 

The operators Z k , Z k with k > 1 also obey the contracted Dolan-Grady relations. Hence, 
one can try to use them to define a "new" physical system with the Hamiltonian of the 
same form ( |2.1| ). Since the operators T> and T> are the tensors of helicity —1 and +1, the 
generating elements Z kl Z k of the "new" system can be treated as a deformation of the 
original £7(l)-connection. As a result, the "new" system constructed in such a way turns out 
to be equivalent to the "old" one up to a redefinition of the parameters Co, c\ and f3. 

Among the central charges (|5.2|), besides the Hamiltonian Ti n = |j7^, the operator J 2 is 
the only independent integral of motion, which can be represented in the form 

J 2 = - Af3 2 (2Jl + W) - \{^ 2 ?S* {S* + na 3 ) - \c. (5.10) 
Hence, in the representation under consideration, the Hamiltonian 7i n and the Casimir 



operator W form the set of independent central charges of the nonlinear superalgebra |L5 
which are the differential operators. We also note that the obvious integrals of motion, S and 
(T3, of the super symmetric system can be transformed into the central charges T\ = 2S + na 3 
and T 2 = S(S + na^), and in terms of them, the second structure in the integrals ( |5.10|) is 
represented linearly, 

5$ (5$ + na 3 ) = - i$Ti + T 2 . 

All the commuting charges J l n (|5~2^ ) with / > 3 are not independent since they can be 
represented as: 

Thus, from the above algebraic analysis it follows that the operators W, T\ and T 2 
exhaust the list of independent nontrivial central charges of the n-HSUSY system described 
by the Hamiltonian 7i n and by the supercharges ( p,12|) . The given systems on the sphere 



and hyperbolic plane provide the first examples of the nontrivial realization of the contracted 



Onsager algebra in contrast with the cases of the plane [13] and ID systems with n-HSUSY 



T2|| , where the realization of the algebra turns out to be trivial [|T^ . 
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In conclusion of this section, let us return to the integral of motion J = W/VC con- 
structed from the Casimir operators of the intrinsic algebra. Its coordinate form is 

J = zd - zd - S*. (5.11) 

Since the J is associated with the axial symmetry of the system, it has a sense of the 
total angular momentum operator. Then, the shifted helicity S<$> ( |5.4T) can be treated as 
the "effective spin" , controlling the phase of the quantum wave functions that they acquire 
under the 27r-rotation of the system. In the spherical case, its eigenvalue s$ (|2.21|) is half- 
integer when the quantized magnetic flux $ takes an odd value, i.e. we have here some 
kind of bose-fermi transmutation. On the other hand, in the case of the hyperbolic plane 
the magnetic flux is not quantized, the s$ can take arbitrary values, and one can say that 
the quantum states of the system are anyonic-like. However, such an interpretation has a 
weak point: the operator J was obtained by us from the intrinsic algebra, whose iso(3) (or 
iso(2, 1) in the hyperbolic case) generators E±, T±, To are not integrals of motion, and, so, do 
not correspond to any symmetry of the system. Therefore, it seems that nothing prohibits 
to shift ( |5.11| ), say, by a constant term /(co,Ci) such that /(0, 0) = 0. Below, however, we 
shall give an additional argument in favour of treating (|5.11|) as a total angular momentum 
operator. 



6 Algebraic approach to analysis of symmetries in an 
external field 

In the previous section we have established that the integral of motion corresponding to 
the axial symmetry of the system is the Casimir operator of the intrinsic algebra generated 
by the covariant derivatives ( |5.3|) . As a consequence, the operator J = W/C automatically 
commutes with the Hamiltonian constructed in terms of covariant derivatives. In contrast 
with the generic case ( |2.15|) , the constant {c\ = 0) magnetic field does not break the isometry 



of the sphere (hyperbolic plane). Therefore, the system defined by the Hamiltonian ( 2.1 
with B = const has to have three integrals of motion corresponding to the isometry of the 
background space. Here we analyse the symmetries of the system in the constant external 
magnetic field in the context of the corresponding intrinsic covariant algebra. This will help 
us to understand better the question of fixing the form of the total angular momentum 



operator (|5.11 ). Moreover, our analysis will result in formulation of the general method 



(beyond the context of supersymmetry) of obtaining the explicitly covariant form of the 
integrals of motion. 



Using the representation (|2.17 ) of the "bosonic" Hamiltonian in the gauge ( 2.18|) with 



Cx = 0, one can find the coordinate form of the three integrals of motion, 

J+= - (3 2 z 2 d -d+\2s-^j (3 2 z, 

J = zd - zd - s + — , (6.1) 
J_ = d + p 2 z 2 d - |/5 2 z, 
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where we assume that these operators act on the tensors of helicity s. Here, the constant 
term in the operator J is completely fixed by the condition that the operators (|6.1| ) realize 
a unitary representation of the non-Abelian algebra su(2) being the isometry algebra of the 
sphere. The modification corresponding to the hyperbolic plane is reproduced, as always, 
by the formal change (3 —>■ i/3, and in this case the integrals span the algebra su(l, 1). The 
detailed discussion of different aspects of the systems on the sphere and hyperbolic plane in 
the constant magnetic field can be found in Refs. |35|, |3~I, 



The integrals ( |6.1| ) are given in the non-covariant form. Nevertheless, we know that 
they have to be the scalar operators since we consider the spectral problem on the tensors 
of the definite rank, s, and the symmetry, generated by the integrals, corresponds to the 
degeneracy of the spectrum. The coordinate form ( |6.1|) and the information about the type 
of the operators is enough to properly define their behaviour on the (co)tangent bundle, 
but the covariance with respect to the corresponding [/(l)-bundle is still vague. Below, 
modifying appropriately the algebraic method used in the previous section, we develop a 
general method of finding the explicitly covariant form of the integrals. 

Let us consider the algebra generated by the operators Z and Z in the case of the constant 
magnetic field, 

[Z, Z] = -4(3 2 S^ } [<S$, Z\ = -Z, [5$, Z] = Z. (6.2) 



Here the operator 5$ is formally the same as for the algebra (|5.3| ), and (3 is assumed to 
be pure imaginary in the case of the hyperbolic plane. One can note that this algebra is 
exactly the isometry algebra of the background space. As we will see later, this is a general 
(local) property of the spaces of constant sectional curvature when external fields are absent 
or constant. 

Unlike the case of the inhomogeneous magnetic field ( |2.15[ ), in the algebra ( |6.2| ) gener- 



ated by the covariant derivatives only one additional element appears. Obviously, we cannot 
construct the integrals ( |6.1D only in terms of the operators Z, Z and S<$>. In the inhomoge- 
neous case the components of the Killing vector corresponding to the axial symmetry were 
the necessary constituents of the covariant construction of the integral W ( |5.5| ). For the 
homogeneous magnetic field one needs to use the three Killing vectors. The form of their 
components in the coordinate system we are working in is the following: 

V = ^p(zz)z, V = z, V^^zz)^ 2 , V 1 = p 2 z 2 , V 2 = ^p(zz), V 2 = 1. 

Under the change of the coordinates the V, {V), i = 0,1, 2, are transformed as the tensors of 
helicity —1 (+1). We will treat Vi and V as multiplication operators for the given quantum 
system. In this context, the operators Vi and Vi are mutually adjoint with respect to the 
scalar product Q2.6|) . Together with the covariant derivatives, Vi and Vi generate the following 
algebra: 



(6.3) 



[Z, Vo] = 


Go, 


[z, 


V ] = 


—Go, 


[Z, Go] 


= -4/3X 


[Z, Go] = 


4{3 2 Vo, 


[Z, V 1 ] = 


G\, 


[z, 


Vi] = 


—Gi, 


[Z, G x ] 


= 4/3V 2 , 


[Z, G 1 ] = 


4P 2 V X , 


[Z, V] = 


-Gi, 


[z, 


v 2 ] = 


Gi, 


[Z, G x ] 


= -4/3V 1; 


[Z, GJ = 


-4f3 2 V 2 , 



[5*, V] = -V u [S , Vi] = V u 
[Z, Vi] = [Z, Vi] = [5*, G a ] = [5*, G a ] = 0, 
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where i = 0, 1,2, a = 0, 1, and in the chosen coordinate system, G = G , G x and G\ are 
the scalar multiplicative operators of the form 

G = ly/p^z) - 1, G x = 2/5 2 z v / p(ii), G x = 2(3 2 z^/^). 

Since Vi, V, G a and G a are the multiplicative operators, they form an Abelian subalgebra. 
Moreover, from the commutation relations ( |6.3| ) it follows that this Abelian subalgebra is 
an ideal. Finally, one can conclude that the algebra given by the commutation relations 
(§3), Q is = (t_(3) © to (3) © t + (3))© s isom, where t_(3) = span^i, G x , V 2 }, t (3) = 
span{Vo, G ,V }, t+(3) = span{V2, G x , V x }, and isom is the isometry algebra of the sphere or 
of the hyperbolic plane. The Abelian subalgebras t±(3) are mutually conjugate with respect 
to the scalar product ( |2.6| ). The algebra g has seven Casimir operators constructed from the 
generators t±(3), to (3). They are trivial in the given representation and their particular form 
is not important for the further discussion. 

The Hamiltonian ( |2.1|) is composed of the covariant derivatives, which belong to the 
algebra isom. Therefore, the Casimir operators of the subalgebras t±(3) ffi s isom and to (3) © s 
isom, each of which is isomorphic to iso(3) (sphere) or to iso(2, 1) (hyperbolic plane), 

J + = ZV 2 -V X Z+(S*-1)G X , 

J = V Q Z + ZV - (S* - 1) G , (6.4) 
J_ = V 2 Z-ZV X + (S*-1)G X , 

commute with the Hamiltonian ( |2.1| ). These scalar operators are explicitly covariant with 
respect to the U(l) and (co)tangent bundles and in the gauge ( |2.18|) their coordinate form 
coincide with the operators flfi.lp . Thus, we have found the covariant form of the integrals 
of motion for the given system (with B = const) in the pure algebraic manner. 

It is worth noting that in the given representation the multiplicative operators obey some 
additional relations of a polynomial form. We will discuss only the most interesting of them. 
For example, one can verify that the multiplicative operators satisfy the constraints 

2/3V - G 1 V 1 - G 1 V 2 = 0, (1 - G )V 2 - G X V = 0, (1 + G Q )V X - G X V = 0, 

and their conjugate. These constraints mean that the Killing vectors are not independent 
in some sense. The commutation of these relations with the covariant derivatives leads to 
another set of constraints: 

G a = A/3 2 (V V 2 + VlVq) , G a = 4(3 2 (VbK + V 2 V ) , G = 2 (V 2 V 2 - V X V X ) . 

Obviously, one can use these relations to exclude the operators G , G x and G x from the 
algebra (|6.2j) , (|6.3| ). Then, one can say that the operators Z, Z, Vi and Vi generate some 
quadratic algebra. Therefore, the quantum system in the homogeneous magnetic field rep- 
resents a physical system with the intrinsic nonlinear algebra. 

In the case of an inhomogeneous magnetic field with the axial symmetry, the components 
of only one of the Killing vectors, say, V and V , can be considered as the proper multiplica- 
tive operators of the corresponding quantum system. Thus, for the axial magnetic field of a 
general form only the operators Z, Z, Vq and Vq can be considered as generating elements of 
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the algebra, intrinsic to the given quantum system. In general case, the covariant derivatives 
Z and Z generate an infinite-dimensional subalgebra. Nevertheless, the class of systems with 
a finite-dimensional generated algebra should be large enough. Note, that the appearance of 
such a finite-dimensional algebra implies that the magnetic field obeys some conditions (like 
those appearing from the Dolan-Grady relations). The algebraic approach for obtaining the 
corresponding integrals of motion can be applied in these cases as well. 

In addition to the gauge interaction, one can include into this scheme the interaction 
with a scalar potential, which, in general, should be treated as an independent generating 
element of the intrinsic algebra. As a toy example, one can include into the Hamiltonian 
( |2.1D the interaction with the scalar potential of the form U = U(G , VoVo). Obviously, such 
a potential does not change the intrinsic algebra g ( |6.2|) , ( |6.3D , but it breaks the symmetry 
of the quantum system to the axial subgroup generated by the operator Jo- 
in the spherical case, when the magnetic field is switched off, the generators (|6.1|) , forming 
the su(2) algebra, realize the representations with integer weight only since we start form 
the wave functions to be tensor fields, i.e. initially spin of the system is integer. When the 
magnetic field is switched on, in accordance with the coordinate form of Jq (|6 . 1|) , the spin 



of the quantum system effectively undergoes the shift by <3>/2. Hence, if the flux parameter 
$ is odd, the spin of the system is half-integer. This phenomenon is well-known, e.g., in the 
context of the charged system immersed into the field of a magnetic monopole {|7|]. At the 



same time, the 2D quantum system on the sphere subjected to the homogeneous magnetic 
field can be treated as a 3D system in the field of magnetic monopole reduced to the spherical 
geometry. In the case of the hyperbolic plane the real parameter $ is not restricted, and 
the generators (|6.1|) forming the noncompact su(l, 1) algebra realize a representation with 



arbitrary real weight [^8J defined by the shifted helicity operator S$. 

The magnetic field of the form (|2.15| ) possesses only the axial symmetry associated with 
the Killing vector Vo, Vq. It means that the components of the other Killing vectors cannot 



be considered as a proper generating elements any more. Moreover, in the case (|2.15|) the 



multiplicative operators Vo, Vo are not independent generating elements since the intrinsic 
algebra ( |5.3|) is generated merely by the covariant derivatives. Indeed, the operators T> and 
T> in the algebra (|5.3j ) are proportional, correspondingly, to Vo and Vq. The integral W ( |5.5|) 
associated with the axial symmetry, being the Casimir of the algebra t(3)© s isom, is similar to 
the integrals ( |6.1| ). In this context, one can say that the presence of the inhomogeneous part 
in the magnetic field ( [2.15|) breaks the algebra g (|6.2|) , (|6.3j ) to its subalgebra to (3) @ s isom, 



and from the three Casimir operators fltt.lp only the operator J survives in the form of the 
integral (|5 . 1 1|) . Such a relation between the cases with inhomogeneous axially symmetric 



magnetic field and homogeneous magnetic field supports our treatment of the operator ( [5.1 1| ) 
as the 2D total angular momentum operator of the system (|2.1|), (|2~5|) , (|2.15| ). 

The above algebraic approach can be generalized to the case of curved spaces of higher 
dimensionality. To argue this, let us discuss a quantum system in a D- dimensional Rieman- 



nian space with a constant sectional curvature ||49|| , the wave functions of which are covariant 



tensors. Then the Riemann tensor has the form 

where the constant 1Z is the scalar curvature. Obviously, such a space is a symmetric one 
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and its isometry algebra is maximal: so(D + l) (TZ > 0), iso(D) (71 = 0) or so(D, 1) (TZ < 0). 
The covariant derivatives with the Riemann connection generate the corresponding isometry 
algebra. Indeed, their commutator can be represented in the form 

[V„ V,] = l -R^M aP = K—Mp,, 

where the operators M^ v are given in a matrix representation defined by the fields which 
the covariant derivatives act on. For example, the action of the operators on a contravariant 
vector field reads as o V a = (5®g uy — 5£<7 M7 )V 7 , while for a spinorial field the corre- 
sponding matrix is proportional to [7^, j v ], where 7^ are the D-dimensional 7-matrices. The 
generalization to the fields of arbitrary type is straightforward. It is not difficult to verify 
that the covariant derivatives and the operators form the algebra isomorphic to the 
isometry algebra of the background space, i.e. the covariant derivatives are the generating 
elements of the isometry algebra. The number of Killing vector fields V a = V^(x)d^ forming 
the isometry algebra is equal to D(D + l)/2. Let us suppose that the interaction with an 
external gauge field or with a scalar potential breaks this algebra to a subalgebra spanned by 
a subset of Killing vectors Vi, i — 1, . . . , I < D(D + l)/2. Treating the components V? and 
the scalar potential (if it is present) as multiplicative operators of the given quantum system, 
one can consider them together with the covariant derivatives as generating elements of the 
intrinsic algebra of the system. The corresponding Hamiltonian is supposed to be composed 
of the covariant derivatives and of the scalar potential. Evidently, if the resulting algebra is 
finite-dimensional, then the covariant form of the corresponding integrals of motion can be 
found in closed terms. 



7 Discussion and outlook 

Let us summarize briefly the obtained results and discuss some open problems that deserve 
further attention. 

We have discussed the application of the general algebraic scheme of the nonlinear holo- 
morphic supersymmetry to the quantum mechanical systems with nontrivial 2D Riemann 
geometry in the presence of external magnetic field. The analysis of Dolan-Grady relations 



[2^, underlying the construction [15|, shows that 



• The nonlinear holomorphic supersymmetry can be realized on Riemann surfaces with 
constant curvature only. 

Moreover, for non- vanishing curvature the n-HSUSY can be realized solely in the case of the 
contracted Dolan-Grady relations (|2.11 ). We have investigated the cases of the sphere and 



Lobachevski plane, which together with the zero curvature case of the plane analysed by us 
earlier |L3|], exhaust the simply connected Riemann surfaces. We demonstrated that for the 
n-HSUSY system on the sphere and hyperbolic plane a part of the spectrum can be found 
algebraically by reducing the corresponding 2D Hamiltonian to one- dimensional families of 
the QES Hamiltonians described by the s[(2,M) scheme. 

We have found that the one-dimensional QES potential corresponding to the sphere 
admits a discrete group of transformations of the parameters of the potential. The ID 
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potential corresponding to the case of hyperbolic plane admits a similar discrete group of 
transformations if some parameter is restricted to be integer. In particular, this means that 



The discrete transformations fl4.1| )-( fO| ) partially relate distinct sectors of the spectrum 



of the same 2D system, or the spectra of different 2D systems. 

Due to the second relation in ( |4.4j ), the discrete transformations are a kind of duality trans- 
formations. We argued that, in general, other QES sl(2, K) potentials also can admit such 
a symmetry group after discretization of some parameters. This was illustrated by the ex- 
ample of the QES potential (|3.9| ). We are going to investigate in detail the revealed discrete 
symmetry of QES potentials elsewhere. 

The analysis of the algebraic content of the model on the sphere and hyperbolic plane 
shows that all the observables of the systems can be represented in terms of the generators of 
the iso(3) (sphere) or tso(2, 1) (hyperbolic plane) intrinsic algebras. All the higher generators 
of the contracted Onsager algebra associated with the relations Q2.ll ) are linearly dependent 
(see Eq. (|5.9|)) but non- vanishing. Hence, 



The supersymmetric 2D systems on the sphere and hyperbolic plane realize the non- 
trivial representations of the contracted Onsager algebra. 



Using the representation (|5.9|) , we have found that besides the Hamiltonian, the infinite set 
of commuting charges (|5.2|) contains, in fact, only one non-trivial independent integral of 
motion, which is a differential operator being, simultaneously, the Pauli-Lubanski "pseudo- 
scalar" W (|5.5|) of the corresponding i5o(3) or iso(2, 1) algebra. It plays the role of the 
2D total angular momentum ( |5.11| ), and its noncoordinate part has a sense of the effective 



spin of the system. The value of the effective spin 5$ (|5.4j ), being a helicity shifted by the 
magnetic flux parameter, takes integer or half-integer values on the sphere, and corresponds 
to the anyonic-like quantum states on the hyperbolic plane. 

The algebraic analysis allowed us to find the covariant form of the integral of motion 
associated with the axial symmetry. Following the same pattern we have found the covariant 
form of all the integrals in the case of the homogeneous magnetic field. On the basis of these 
explicit examples, we have proposed 

• A general algebraic method to find the covariant form of integrals of motion of a 
quantum system in a curved space in the presence of an external gauge field. 

The integrals obtained in this way have transparent commutation properties. We hope that 
this approach will be helpful for a wide class of quantum systems with symmetry. 

Nowadays, a considerable attention is paid to a generalization of the 2D systems on Rie- 
mann surfaces with a constant magnetic field |3l], [32], |33|, |35| to the case of noncommutative 



surfaces |50], [51], 0, Such quantum systems are exactly solvable in the both cases of 
commutative and noncommutative spaces. The supersymmetric systems considered in this 
paper can be treated as a generalization of the constant magnetic field case to the case of 
the nontrivial field (B ^ const) accompanied by the transition from exactly solvable systems 
to quasi-exactly solvable ones. Therefore, the interesting problem is to generalize the 2D re- 
alizations of the nonlinear holomorphic supersymmetry to the cases of the noncommutative 
plane, sphere and hyperbolic plane. Such a generalization will be presented elsewhere [53] . 
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